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■ ABSTRACT 

m ■ 

, The spin-down power of a pulsar is thought to be carried away in an MHD wind in 

,—1 \ which, at least close to the star, the energy transport is dominated by Poynting flux. 

\ The pulsar drives a low-frequency wave in this wind, consisting of stripes of toroidal 

i-^^' magnetic field of alternating polarity, propagating in a region around the equatorial 

I ' plane. The current implied by this configuration falls off more slowly with radius than 

2 ' the number of charged particles available to carry it, so that the MHD picture must, at 

c/3 . some point, fail. Recently, magnetic reconnection in such a structure has been shown to 

• • , accelerate the wind significantly. This reduces the magnetic field in the comoving frame 

• ^ i and, consequently, the required current, enabling the solution to extend to much larger 

^ I radius. This scenario is discussed and, for the Crab Nebula, the range of validity of the 

' MHD solution is compared with the radius at which the flow appears to terminate. For 

sufficiently high particle densities, it is shown that a low frequency entropy wave can 
propagate out to the termination point. In this case, the 'termination shock' itself must 
be responsible for dissipating the wave. 
Dedicated to Don Melrose on his 60th birthday 
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1. Introduction 



The supply of relativistic electrons and magnetic field needed to provide the synchrotron 
radiation observed from the Crab Nebula has long been thought to originate in a central star 
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(Piddington 1957; Kardashev 1964). With the identification of this star as a pulsar, the earher 

suggestion (Pacini 1967) that the star also emits a wave, whose energy is presumably released 
into the Nebula, attracted considerable attention (Ostriker &: Gunn 1969; Rees &: Gunn 1974). 
As an energy source for its environment, the Crab pulsar is by no means unique. The pulses of 
electromagnetic radiation emitted by other pulsars contain only a fraction of the power which the 
neutron star is apparently releasing from its store of rotational energy. In several cases, nebular 
emission is observed to surround the pulsar, but, even where it is not, it is generally thought that 
most of the spin-down power is carried away from the star by a relativistic wind consisting of a 
mixture of particles, waves and magnetic field. Attempts to provide a consistent description of this 
wind touch upon a fundamental problem in electrodynamics — the dichotomy between a single 
particle approach and a continuum description. In the pulsar wind case, one can either start from 
the solution describing a rotating magnetised object, compute individual particle trajectories and 
treat their combined effect as a perturbation, or, alternatively, consider the particles as a fluid and 
construct, for example, a solution with an ideal MHD wind emerging from a magnetised rotating 
star. 

In the former case, it has been found that the wave component of the vacuum fields transfer 
energy rapidly to the particles. The resulting damping is strong and the wave will not propagate 
if the particle density exceeds a critical value (Asseo et al 1978; Leboeuf et al 1982; Melatos & 
Melrose 1996). Since the particle density decreases as the radius increases, the vacuum wave-like 
solution exists only outside a critical radius. In the fluid case, one faces the difficulty of formulating 
a generalised Ohm's law. The simplest choice of infinite conductivity corresponds to ideal MHD, 
but has the problem that it may imply currents too large to be carried by the finite number of 
charged particles in the fluid. Given the structure of the magnetic field, the current densities can be 
computed from Ampere's law. These generally decrease with radius less rapidly than the decrease 
in the density of charged particles available to carry the current (Usov 1975; Michel 1982). Since 
the velocity of the current carriers cannot exceed that of light, the MHD approach is valid only 
inside a critical radius. Thus, there are theoretical methods to describe both the inner and the 
outer parts of a pulsar driven wind. However, as yet, no way has been found of linking these two 
regimes. 

In the case of the Crab Nebula, observations suggest that the wind energy is dissipated into 
relativistic particles at a radius r ~ lO^rL, where tl = c/Q is the radius of the light cylinder and 17 
the angular velocity of the neutron star. The physics of this region (which we will term 'termination 
shock') is clearly an important problem. As a preliminary step, it is of interest to examine whether 
the wind solutions enable one to make a statement about the properties of the flow immediately 
before entering the shock region. This is the problem we address in this paper. We discuss a highly 
simplified non-axisymmetric ideal MHD wind — a 'striped wind' — and describe how a model 
of magnetic rcconncction (Coroniti 1990; Michel 1994) modifies the wind dynamics (Lyubarsky & 
Kirk 2001, henceforth LK). As a result, the range of validity of the MHD solution is extended. 
Applying the results to the Crab pulsar wind, we find that for sufficiently high particle densities 
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the solution is valid up to the location of the termination shock. We then speculate that the crucial 
physical process at work in the shock is the dissipation of the wave energy carried by the MHD 
wind. 

2. The MHD wind 

An exact analytic solution for an MHD pulsar wind has been found only for the idealised case 
of a monopole magnetic field in the force- free limit (Michel 1973). This solution has a very simple 
structure: the flow does not collimate, i.e., both the velocity and the magnetic field have purely 
radial components in the meridional plane, Bg = 0, = and there exist no closed field lines. A 
more realistic case would allow for closed field lines near the star, but, outside this region, it might 
be expected that the flow along open flcld lines mimics the monopole case, perhaps with a modified 
dependence on latitude. The obvious inconsistency of using a magnetic monopole as a source can be 
lifted by introducing the 'split monopole'. In this case the system retains axial symmetry; the sign 
of each magnetic fleld component changes at the equator of the star, all other quantities remaining 
the same. A current sheet is then implied, which lies in the equatorial plane and separates the 
oppositely directed toroidal field components in the northern and southern hemispheres. It is also 
possible to lift the force-free approximation, replacing it by an ultra-relativistic approximation 
(Bogovalov 1997) to find that highly super- Alfvcnic, rclativistic winds also collimate only very 
slightly and the radial velocity pattern remains a good approximation, even when particle inertia 
is taken into account (Beskin et al 1998; Chiueh et al 1998; Bogovalov & Tsinganos 1999). 

In a radial wind, the outward pressure gradient exerted by the magnetic field is precisely 
compensated by the inward tension force, so that the wind speed remains constant, as does the 
quantity a, defined as the ratio of the Poynting flux to the energy flux in particles. Thus, a 
steady MHD wind driven by a magnetised neutron star carries energy predominantly in the form 
of Poynting flux, which is not converted into kinetic energy anywhere in the flow. This is a major 
difficulty, because modelling of the synchrotron emission of the Crab Nebula (Rees &: Gunn 1974; 
Kennel & Coroniti 1984; Emmering &; Chevalier 1987) suggests that if the wind terminates at a 
rclativistic MHD shock, the energy flux must be carried predominantly by particles. 

The rotating neutron star of a pulsar does not have an axisymmetric magnetic fleld, so that 
an axisymmetric wind is quite possibly a poor approximation. Non-axisymmetric solutions are, 
of course, very much more difficult to flnd. Fortunately, however, the technique of replacing a 
monopole source by a split monopole can also be Tised when the magnetic axis does not lie along 
the rotation axis (Bogovalov 1999). The solution is essentially the same as in the monopole case: 
the velocity field is unchanged, the pattern of the magnetic field lines is the same, but the direction 
(sign) of the field depends on the hemisphere in which the field line has its anchor point on the 
star. The current sheet, which lies in the equatorial plane in the axisymmetric case, takes on the 
form of an outward moving, spiral corrugation, as shown in Fig. 1. 
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Fig. 1. — The current sheet in the spht-monopole solution for a radial relativistic wind. The sheet 
has been truncated at a distance from the axis of 40rL, and a sector has been removed for clarity. 
The wavelength of the corrugation is 2iTri^. The upper figure displays the sheet for an inclination 
angle x between the rotation and magnetic axes of 6°, the lower for an inclination angle of 26°. 
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In the ultra-relativistic, highly super- Alfvenic Umit, the solution in terms of spherical polar 
coordinates (r, 9, (p) is given by: 

Ur = /^TL U0 = O U<p = 

Br = ±BQ{n^/rf Bg = B^ = ^{rL sin e/pr)Bo (1) 

^^ = E0 = TBosme{rL/r) = 

(Bogovalov 1999), where «r,6»,(/) are components of the four velocity, B and E are the magnetic and 
electric fields and Bq is the magnitude of the radial component of B at the light cylinder, r = tl- 
The signs of Br, B^ and Eg depend on whether the field line lies above (upper sign) or below (lower 
sign) the current sheet. Associated with this solution are distributed ('body') currents and charges 
as well as surface currents and charges located in the sheet. The body current is found from the 
curl of B, (since the displacement current vanishes outside the sheet) and the body charge follows 
from Poisson's equation: 

> = -(c5oCOS0/27r/3rL)(rL/r)2 je = = 

P = f3jr/c ^ ' 

The current sheet is located at the position $(r) = 0, with 

/ r — j3ct \ 

$ = cos X cos ^ + sin X sin 9 cos ( — ^ h (pj (3) 

where x is the angle between the magnetic symmetry axis and the rotation axis and /3c ~ c is the 
3- velocity of the wind. The surface current carried in the sheet is directed perpendicular to the 
adjacent magnetic field, and has the magnitude 

, c[l-/32(n.f)2] ^, 

k = — ^ |n X B| (4) 

where n = V$/|V$| is a unit vector normal to the sheet [e.g., Jackson (1975), pp 20-22]. Except 
at points where the sheet normal is perpendicular to f , the surface current scales with 1/r, whereas 
the body current scales according to jr oc 1/r^. The particle density in the wind drops off as 
1/r^, so that the charge carriers can maintain the body currents without an increase in velocity. 
To maintain the surface currents, however, an increase with radius of either the surface density of 
charge carriers, or of their velocity is required. The MHD solution fails if equation (4) requires a 
current too large to be carried even if the particles contained in the sheet move at the speed of 
light. 



3. Reconnection 

At large radius, the corrugated current sheet described above resembles locally, a set of con- 
centric nested spheres, when viewed at latitudes close to the equator (|90° — 0| <^ %). The dynamics 
of the sheet may then be considered using a one-dimensional (radial) treatment. In the equatorial 



- 6 - 




2n in 

r/ri^ — 27vN 

Fig. 2. — A schematic sketch of the radial structure of the magnetic field and current density in the 
wind far outside the light-cylinder (A^ ^1). In the equatorial plane, the current sheets are equally 
spaced, as shown. However, except in the singular case x = 90°! the intersheet spacing depends on 
the magnetic polarity away from the equator and the sheets vanish altogether for |90° — 0\ > x 
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plane, the sheets of opposite polarity are evenly spaced. A sketch of the configuration at r » is 
shown in Fig. 2. 

According to Ampere's law, the integral of the current density across a current sheet is propor- 
tional to the change in magnetic field strength, measured in the comoving frame. In a radial wind, 
the toroidal field drops as 1/r, and the particle density as when both quantities are measured 

in the lab. frame. Thus, in a wind of constant speed, the charge carriers in the sheet are forced 
to counter-stream with higher and higher velocities as the plasma moves outwards. This situation 
is likely to lead to instability, which could, in turn, provide an anomalous resistivity and hence 
magnetic reconnection. Coroniti (1990) [see also Michel (1994)] proposed a simple picture intended 
to capture the physics of magnetic reconnection in the sheets. The suggestion is that the sheets 
attain a minimum thickness equal to the gyro-radius which a particle in the hottest part of the 
sheet would have, if it were moving in the magnetic field adjacent to the sheets. Equivalently (to 
within a factor of order unity), the thickness can be set by demanding that the maximum current 
density in the sheets is equal to the product of the particle density, the electronic charge and the 
velocity of light. This enables one to formulate a system of equations for a wind that consists of 
two phases: a hot unmagnetised phase (corresponding to the current sheets) and a cold magnetised 
phase. From Fig. 1, it can be seen that, for an oblique rotator, the points at which a radius vector 
cuts the current sheet are equally spaced only in the equatorial plane. At other latitudes, the cold 
phase of one polarity dominates. In order to maintain pressure equilibrium, the magnetic field 
strength in the cold phase must remain constant on the scale of a wavelength, so that the cold 
phase of the dominant polarity is thicker than that of the opposite polarity. At a critical latitude 
{9 = 90° — x) the radius vector grazes the current sheet, and for |90° — d\ > X the hot phase is 
absent. 

In Coroniti's original work, the density in the hot and cold phases was assumed equal. How- 
ever, this leads to inconsistencies and an incorrect picture of the evolution of the wind. Recently, 
Lyubarsky & Kirk (LK) have used a short wavelength approximation to analyse the evolution of the 
pattern of hot and cold phases, which is simply an entropy wave in the wind plasma. Reconnection 
at the sheet boundaries starts at a certain radius r = rgtartj which depends on how many particles 
are contained in the current sheet initially. The short wavelength approximation consists in assum- 
ing this radius to lie well outside the light-cylinder: e = TL/rstart ^ 1- The slow evolution of the 
wave is governed by a set of five equations. Three of these are conveniently written in differential 
form. These are the conservation of particle number and of energy. 




(5) 



(6) 



and the entropy equation 



B? dR 
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Here, 70 and vo{R) are the zeroth order (in terms of e) Lorcntz factor and (3-) velocity of the 
flow, the dimensionless radius variable is i? = r/rstart) ^(R) is the fraction of a wavelength occupied 
by the sheets, n'^{R) and are the (zeroth order) proper number densities in the cold and hot 

phases, respectively, and B{R) = ^qB'{R) is the (zeroth order) magnetic field in the cold phase 
(the prime indicates a quantity measured in the plasma rest frame). The condition of ideal MHD 
(flux freezing) in the cold phase: 

B 

= constant (8) 

rric 

follows from the zeroth order equation of continuity, r'^v^n'^ = constant and the zeroth order Faraday 
equation, \^rE = —^-^B = together with the ideal MHD condition E = —vB. The remaining 
equation stems from the requirement that the sheet thickness equal the gyro radius: 

rn'c , , 

^ = , / 9 

In this equation, we have introduced the multiplicity parameter k, which is a dimensionless measure 
of the number density in the cold phase of the wind, in terms of the Goldreich-Julian density (see 
LK): 

' = B^ 

where P is the pulsar period, ni, = rii/e^ and 5l = Bi/e are the density (in the cold phase) and 
magnetic field, extrapolated back to the light cylinder, in terms of the values ni, Bi aX r = rgtart- 

The set of differential equations (5), (6) and (7) together with (8) and (9) have been solved 
by LK. The solution is specified by three parameters, for example, the Lorentz factor 7l, the 
multiplicity k, and the ratio of the particle gyro frequency at r = tl to the rotation frequency 
of the pulsar ui^/Vl. The values of the first two parameters are uncertain, but i^i,/Vt is more or 
less directly accessible from observation and is, to within a factor of a few, equal to the potential 
difference generated across open field lines, in units of the electron rest mass. (For the Crab pulsar 
u}j,/Vt 10^^.) As an alternative to the multiplicity, one may use the magnetisation parameter cjl, 
defined as the ratio at the light-cylinder of Poynting flux to particle energy flux in the cold phase 
of the wind, and related to the above parameters by ctl = ui,/{2K'yi,Q). In addition, the initial 
condition rrstart must be given. 

In Fig. 3 we display an example of a solution appropriate for a pulsar that generates a maximum 
potential equal to that of the Crab pulsar: uJi^/Q.= 10^^. The multiplicity is chosen to be k = 10^, 
and the initial Lorcntz factor is 7l = 100, (corresponding to (Tl = 5000). Reconnection is assumed 
to commence at rgtart = lOOrL, i.e., e = 0.01. As the magnetic flux is dissipated, the hot phase 
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Fig. 3. — A numerical solution of the reconnecting striped wind problem, using parameters appro- 
priate for the Crab pulsar: iOi,/Q, = 10^^, k = 10^, 7l = 100, corresponding to (Tl = 5000, and 
assuming reconnection starts at rgtart = lOOrL- The dashed line shows the asymptotic solution of 
LK. 
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expands and performs work on the wind, causing it to accelerate. The system quickly approaches 
the asymptotic solution given by LK and shown as dashed lines. This solution reads: 



T' 



7 = (11) 



^ = xh^z:— (12) 



ISvrrLWL 



(l 

Stt ) V rL 



where p' = B'^/Stt is the pressure of the hot phase and T' its temperature, in energy units. 



(13) 



4. Validity of the solution 

The ideal MHD part of the solution, in which the radial velocity is constant and no magnetic 
flux is dissipated, is valid, by assumption, for r < rgtart- This range is determined by the surface 
density of charges assumed to be present in the quasi-spherical current sheets when they are first 
established. For r < rstart the current sheet contains a sufficient number of particles to supply the 
required current. Provided the initial surface density is small compared to the surface density of 
particles in the cold, magnetised phase (i.e., the integral of over one wavelength) the sheets are 
able to maintain the required current at r > rgtart by expanding and absorbing particles from the 
cold phase. In this way, the MHD picture of the wind retains validity, although the flow necessarily 
contains non-ideal regions in which oppositely directed magnetic flux is steadily annihilated. The 
solutions can be followed until the current sheets — as seen in the one-dimensional description — 
appear to merge. In the equatorial plane, this occurs at a radius given approximately by 

^max ~ Y^^L (15) 

and implies the complete annihilation of the magnetic field. Above and below the equatorial plane, 
the current sheets arc not equally spaced, and the total magnetic flux passing through one complete 
wavelength of the pattern is non-zero. Since this flux is conserved in our picture, reconnection leaves 
behind a residual magnetic field, whose magnitude can be found from Eq. (3) and is 

_ tl-Bl /2arccos(cotxcot^) 

-Ores ~ I 

r \ TT 

Beyond rmax^ the 'current sheet' is no longer corrugated, but lies in the equatorial plane. Since 
the toroidal magnetic field goes smoothly to zero as this plane is approached, the current it carries 
is proportional to the radial component of the field, which varies as r~^. As in the axisymmetric, 
split-monopole case, the MHD picture is still valid, since the number of available particles is always 
sufficient to provide this current. 
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However, there are two limitations which restrict the vahdity of this solution at large radius. 

The first is the assumption, implicit in Eq. (9), that the electrons in the hot phase arc relativistic. 
This is a relatively minor technical limitation of the treatment, which is encountered for radii more 
than a few percent of rmax- The second limitation is more subtle, and can be stated in a number of 
ways. It also arises from Eq. (9). In prescribing the thickness of the sheet in terms of the densities 
in the two phases and the wind velocity, the movement of the sheet edge is decoupled from the 
local wave speeds, and can even move supcrhiminally. Clearly, a physically realistic prescription 
would involve, at the very least, an additional dissipation of energy into waves when the velocity 
of the sheet edge approaches c. This point is reached when 

Wedge = Vo(l + ) = C (17) 



which can be written, using Eq. (12) 



2r 

A = (18) 

To within a factor of the order of unity, this is coincident with the limit discussed in LK [Eq. (36)]. 
Using the asymptotic solution, the critical distance Vct beyond which the solution fails is 

rcr = 1.5x103kV(^)~' (19) 
Combining this with equation (15), we find 

^ = 1.5xl0V(^^)"' (20) 



5. Discussion 

The inclusion of non-ideal MHD effects — albeit phcnomenologically — results in a substantial 
extension beyond rgtart of the range of validity of the 'inner' wind solution, which can be described 
using a continuum picture. Starting with an entropy wave outside the light cylinder, this solution 
can exist up to a maximum radius ~ rmax (Eq. 15) and convert a substantial part of the total energy 
flux into bulk kinetic energy of the particles, provided k > O.I^/cjl/^^ [see eq. (20)]. The situation 
in the outer parts of the wind then resembles the axisymmctric, split monopole solution, except 
that the magnitude of the toroidal component of the magnetic field passes smoothly through zero 
on crossing the equatorial plane. However, if K<0.1ytUL/O, the solution loses validity at r ~ Vcr 
(Eq. 19) and the energy contained in the non-axisymmetric or (entropy) wave component must be 
channelled into another type of wave and/or into the particles. 

In reality, the wind driven by a pulsar is confined by the surrounding medium, which may 
disrupt the solution either before reconnection is complete or before the MHD approach loses 
validity. Usually, disruption is envisaged to occur at a 'termination shock' where the flow makes 
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an abrupt transition from supersonic to subsonic velocity. In the case of the Crab Nebula, this 
transition can be identified with a sharp increase in the synchrotron cmissivity at a radius of 
approximately ^shock = lO^rL- The Crab pulsar has cjl/^ = 10^^, so that, according to Eq. (15), 
^max. ^ ^shock and Very little of the striped magnetic flux has time to reconnect. For the MHD 
solution to remain valid up to rghock a minimum number of particles is required. Prom Eq. (19) 
this may be expressed as a condition on the multiplicity parameter: k> 1.6 x 10^. Assuming that 
this condition is satisfied, the region r < rghock can be described by the striped MHD wind. Close 
to the equator, most of the energy flux is carried by the magnetic field of this entropy wave. If, in 
this region, the termination shock is to dissipate a substantial fraction of the wind energy, it must 
not only decelerate the wind to subsonic velocity, but also damp the oscillating component of the 
magnetic field. This conclusion is consistent with the constraints from modelling the synchrotron 
emission of the nebula, which indicate that the energy density of the plasma downstream of the 
shock is not dominated by the magnetic field. The physics of such a transition is at present unclear, 
but is likely to have profound implications for particle acceleration. 

We are grateful to L. Ball, Y. Gallant, J. Kuijpers, A. Melatos and D. Melrose for helpful 
discussions. 
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